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Kaon B-parameter using Overlap Fermions 

Thomas DeGrand (MILC collaboration) 

Department of Physics, University of Colorado 
Boulder, CO 80309-390, USA 

I present first results from an in-progress calculation of Bk in quenched approximation using overlap fermions. 
My particular implementation of the overlap uses a kernel with nearest and next-nearest neighbor interactions 
and HYP-blocked gauge connections. Matching to the continuum NDR regularization is done perturbatively. I 
present preliminary results at j3 = 5.9 and 6.1 (lattice spacings 0.125 and 0.09 fm) for quark masses, pseudoscalar 
decay constants, and B-parameter - B^ DR \fi = 2 GeV) ~ 0.66(3 — 4). 



The kaon B-parameter Bk, defined as 
l{m K f K ) 2 B K = (K\s^(l- 75 )ds^{l- l5 )d\K), 
has been computed many times with lattice meth- 
ods. Lattice calculations of Bk require ac- 
tions with good chiral properties, to prevent 
operator mixing with wrong-chirality operators 
from contaminating the signal. There has been 
a continuous cycle of lattice calculations us- 
ing fermions with ever better chiral properties. 
This calculation is yet another incremental up- 
grade, to the use of a lattice action with exact 
SU(Nf) ® SU(Nf) chiral symmetry, an overlap 
action. These actions have operator mixing iden- 
tical to that of continuum-regulated QCD. 

The overlap action used in these studies (!]] is 
built from a kernel action with nearest and next- 
nearest neighbor couplings, and HYP-blocked 
links HYP links fatten the gauge links without 
extending gauge-field-fermion couplings beyond 
a single hypercube. This improves the kernel's 
chiral properties without compromising locality. 
The kernel action is designed to resemble the ex- 
act overlap well enough that its eigenvectors are 
good "seeds" for a calculation of eigenvectors of 
the exact action, and it is kept simple enough that 
finding its own eigenvectors is inexpensive. These 
eigenvectors are used to precondition the calcula- 
tion of quark propagators, in principle eliminat- 
ing all critical slowing down at small quark mass. 

The data set is generated in the quenched ap- 
proximation using the Wilson gauge action at 
couplings (3 = 5.9 (on a 12 3 x 36 site lattice) and 
(3 = 6.1 (on a 16 3 x 48 site lattice) with 40 lattices 



each (so far). The nominal lattice spacings are 
a = 0.125 fm and 0.090 fm from the measured rho 
mass. Propagators for six quark masses are con- 
structed corresponding to pseudoscalar-to-vector 
meson mass ratios of mpsjray ~ 0.6 to 0.85. 



Table 1 

Results from these simulations. 



(3 = 5.9 f3 = 6.1 



1/a (MeV) 


1580(60) 


2190(140) 


TYlnonstrange (MeV) 


4.3(3) 


4.5(3) 


m s trange (MeV) 


105(5) 


110(7) 


m s /m ns 


24.40(4) 


24.41(5) 


U (MeV) 


142(11) 


131(12) 


fx (MeV) 


155(10) 


147(11) 


b (ndr) ^ _ 2 GeV) 


0.66(3) 


0.66(4) 


p>(RGi) 


0.92(4) 


0.93(6) 



The methodology for Bk is well-developed: 
compute an un-amputated correlator which con- 
tains the desired matrix element (two kaon 
sources far apart on the lattice with the four- 
fcrmion operator sandwiched in between), clip off 
the (tukIk) 2 prefactor by simultaneously mea- 
suring the matrix element (0\s~/oj5d\K) , extrap- 
olate/interpolate the lattice B— parameter to its 
value at the kaon mass, and convert the lattice 
number to its continuum-regulated counterpart. 

To maximize the signal volume I computed 
propagators from two well-separated sources 
(N t /2 — 2 temporal sites apart) and brought them 
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together to the operator. I used Gaussian sources 
to maximize overlap onto the ground state. These 
sources do not make momentum eigenstates, and 
so the p = Bk signal is contaminated by a 
p^O contribution. This causes problems at big- 
ger quark mass, because E(p) — m gets smaller 
as the pseudoscalar mass m grows. Fortunately, 
there are two inequivalent paths on the torus to 
disentangle the two "signals," and one can fit the 
Bk correlator to a sum of a p = term and a 
p = 2ir/N s term. This is not a problem at small 
quark mass. 

I extracted a signal from fits to the traditional 
ratio of the K — K amplitude and product of two- 
point functions, as well as correlated fits to the 
K — K amplitude and the two-point functions. A 
ratio plot is shown in Fig. |l|. A typical correlated 
fit is shown in Fig. g, and a plot of lattice Bk 
vs. quark mass is shown in Fig. ||. Uncorrelated 
jackknifc "ratio fits" have small uncertainties and 
are quite stable over a wide range of timeslices. 
However, one would really like to do correlated 
fits. The data points are strongly correlated, and 
it is necessary to use singular value decomposition 
to invert the correlation matrix. When I do this, 
I find that my fits are consistent with the jack- 
knife fits and have reasonable confidence levels. 
To get to the physical kaon I linearly extrapo- 
lated my results with a jacknife; there is no sign 
of discernable curvature in my data. 

I calculated the renormalization factors be- 
tween the lattice- and continuum-regulated 
(NDR) matrix elements using one loop perturba- 
tion theory. At a (lattice spacing) x fi (continuum 
regularization point) = 1, a conversion factor is 
Z = 1 + (a s (q*)/An)z. As one might expect from 
related workQ, the HYP link pushes the constant 
z close to zero. The cost is that q* (defined a la 
Lepage-Mackenzie Q) can move to a small value, 
but sensible values of q* are given by the higher- 
order prescription of || . 

For Bk , the operator (9 + has a matching factor 
Z+ and for the overlap action used here, its pa- 
rameter z = -4.0 at q*a = 0.92 for NDR. The en- 
tire conversion factor for Bk from lattice (3 = 5.9 
to fi = 2 GeV is Z+/Z 2 A = 0.99. (Wilson-action 
kernel overlap actions have z's which are an order 
of magnitude larger.) 
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Figure 1. A traditional ratio plot of the 
Bk graph divided by the product of two point 
graphs, from the = 5.9 data set at quark mass 
am q = 0.050 with axial current sources and sinks 
(squares) and pseudoscalar sources and sinks (oc- 
tagons) . 



How reliable is this number? I have not checked 
it directly (yet), but perturbation theory for the 
local axial vector current can be tested with over- 
lap actions by a comparison of the vacuum-to- 
pseudoscalar meson matrix elements of the axial 
vector and pseudoscalar density. At j3 = 5.9 I 
find Za = 0.97 or 0.98 (depending on the choice 
of lowest-order or higher-order q* , and 0.97(1) 
nonperturbatively. Perturbation theory for the 
matching coefficient fo the MS quark mass can 
also be compared to the nonperturbative calcu- 
lation of Ref. This analysis gives Z(/i = 
2 GeV, a) = 1.10(3) at r Q m PS = 5 and 1.14(11) 
at romps = 3, as compared to the perturbative 
prediction of 0.95. 

My PRELIMINARY values for B K and for 
other relevant parameters (lattice spacings from 
rho mass, decay constants, MS(fi — 2 GeV) 
quark masses) are recorded in the Table. A com- 
bined error from fitting, extrapolation, and lat- 
tice spacing (dominated by statistics) is shown. 



3 




0.9 



0.7 



0.6 



0.5 



i 



[] 

P ~ 



0.0 



0.1 



0.2 



Figure 2. The two axial current correlators (la- 
beled "1" and "2") and the "figure-eight" corre- 
lator (labeled "3"), for the am q = 0.100 ft = 5.9 
data set with axial current sources. A correlated 
fit to three correlators over the range t = 7 — 9 
and 24-28 is also shown. 



My Bk result is in reasonably good agreement 
with the staggered JLQCD result Q a bit higher 
than the CP-PACSQ domain wall fermion re- 
sult and quite a bit higher than the RBCQ do- 
main wall fermion result. It is also consistent 
with the Wilson-overlap results presented by Lel- 
louchjl0| at this meeting. Of course, it is a lin- 
ear extrapolation: that is a dangerous thing to 
do. If the allocation gods allow, I hope to push 
to smaller quark masses (bracketing the neces- 
sary quark mass and possibly revealing the chiral 
logarithm) and collect more statistics (hopefully 
giving more respectable error bars). 

This work was supported by the US Depart- 
ment of Energy. I are grateful to S. Sharpe for 
suggesting this project and C. Bernard and T. 
Blum for helpful instruction. Simulations were 
performed on the Platinum cluster at NCSA. 
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